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ABSTRACT 

The long standing problem of the ordering ambiguity in the definition of the Hamilton 
operator for a point particle in curved space is naturally resolved by using the powerful 
geometric calculus based on Clifford Algebra. The momentum operator is defined to be 
the vector derivative (the gradient) multiplied by — i\ it can be expanded in terms of basis 
vectors 7 M as p = —i^d^. The product of two such operators is unambiguous, and such 
is the Hamiltonian which is just the D'Alambert operator in curved space; the curvature 
scalar term is not present in the Hamiltonian if we confine our consideration to scalar 
wave functions only. It is also shown that p is Hermitian and self-adjoint operator: the 
presence of the basis vectors 7 M compensates the presence of ^J\g~\ in the matrix elements 
and in the scalar product. The expectation value of such operator follows the classical 
geodetic line. 
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1 Introduction 



An important step towards quantum gravity is formulation of quantum field theory in 
curved spacetime [Tj. A commonly used method is De Witt's |2] _ curved space general- 
ization of the Fock-Schwinger proper time technique [Sj for finding the Feynman Green's 
function. In this connection, as a toy model many researchers [2] ^0] have studied the 
quantization of a non relativistic particle in curved space V n , whose dimension and signa- 
ture can be left unspecified. 

The most straightforward approach has been in using the coordinate-based formulation 
of the quantities such as the momentum p^ and the Hamiltonian which is proportional to 
the quadratic form g^p^Pv, (jl, v = 1, ...,n. When p^ are quantum operators the position 
of the curved space metric g^ u (which depends on x^) does matter: the expression p^g^ u p u 
is different from g^ v p^p v or p^Pvg^- This is the well known ordering ambiguity. Different 
ordering prescriptions lead to different Hamilton operators which in general contain a 
term with the scalar curvature ah 2 R. Different authors, using different procedures, have 
obtained different values for the coefficient a. A related problem is the correct definition 
of momentum operator p^ which should be Hermitian and retain the desirable properties 
under general coordinate transformations of coordinates x^. A resolution of the latter 
problem has been proposed by DeWitt [21 El- But the ordering ambiguity has remained 
unresolved within the coordinate based operator formalism. 

A promising more general line of research is geometric quantization jTTj which employs 
coordinate-free differential geometry. If applied to curved spacetime there is no difficulties 
in defining a self-adjoint momentum operator and no ordering problem does arise. 

In this paper we explore the quantization in curved spacetime by using yet another, 
also very powerful, geometric calculus which is based on Clifford algebra ^2j- We define 
momentum as a vector p = 7^(9^ which is expanded in terms of the Clifford numbers 7^ 
which serve the role of basis vectors. The presence of the basis vectors 7^ besides the com- 
ponents p^ renders the corresponding quantum operator and its powers unambiguously 
defined and Hermitian. There is no ordering ambiguity in this approach. We find that 
after quantization the Hamilton operator (entering the Schrodinger equation for scalar 
wave function) contains no term with the curvature scalar R. This result agrees with 
that by Kleinert [7] who has carefully examined the paths integral quantization in curved 
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spaces in a way which is distinct from other approaches found in the literature [HI |H] . 

The quantization based on geometric calculus opens an interesting possibility to define 
the expectation (p) for momentum operator. We show that (p) is real and that it is tangent 
to a geodetic trajectory in curved space. 

In Sec. 2 we briefly review the geometric calculus based on Clifford algebra and the 
concept of vector derivative. In next section we provide a definition of vector integral, an 
operation which is used in Sec. 5 where we discuss the expectation value of the momentum 
operator. In Sec. 4 we consider quantization in curved space and write the matrix elements 
of the momentum vector operator p = 7 M p M between position eigenstates. We find that 
the latter matrix elements satisfy the condition for Hermiticity without adding an extra 
term which is necessary in the coordinate-based approaches [2 -jTU]- We also show that 
(x\p 2 \(j)) can be calculated by inserting twice the complete set of the position eigenstates 
\x) and that the result is the covariant D'Alambert operator (multiplied by —1) as it 
should be. So we have succeeded in providing a consistent and elegant foundations for 
quantum mechanics in curved space. 

2 Geometric calculus 

Since the seminal work by Hestenes ^2|; together with some very important works such 
as those of refs. the geometric calculus based on Clifford algebra has been receiving 
increasing attention |14j . It is a very powerful language for geometry and physics, enabling 
potentially very important generalizations fSJ El EE QB|- The calculus is based on the 
observation that the non commuting numbers satisfying the Clifford algebra relations 

i l ,-iv = \{i^iu + iui^) = g^ (i) 

can represent basis vectors. An arbitrary vector a is then a linear superposition 

a = (2) 

where the components a M are scalars from the geometric point of view, whilst 7^ are 
vectors. The latter relation is just a more general form of a similar relation discussed in 
textbooks: 

a = aol + aVj (3) 
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which states that a 4- vector can be written in terms of the Pauli matrices 0{. Here we do 
not consider 7 M as being necessarily matrices 2 ; they are just numbers satisfying eq.(JIJ. 
Besides the basis 7 M we can introduce the dual basis 7^ satisfying 

/•fEl( 7 Y + 7Y) = f (4) 

where g^ v is the covariant metric tensor such that g^ a g a v = S^ v , 7^7^ + 7^7 M = 25 M v and 
7 M = g^ v ^ v . We shall consider curved space in which 7^ and g, v depend on position x. 
The vector derivative or gradient is defined according to 

d = Yd, (5) 

where is an operator, that I will call geometric derivative, whose action depends on 
the quantity it acts on. 

Applying the vector derivative d on a scalar we have 

d<\> = 7^ (6) 

where <9^0 = (<9/cte M )0 coincides with the partial derivative of 4>. 
But if we apply it on a vector a we have 

da = 7 ^K7„) = Y(d,a^ u + a u d^ lv ) (7) 

In general 7l/ is not constant; it satisfies the relation fF2\ Ho] 

9„iu = r^ 7a (8) 

where V" is the connection. Similarly, for 7^ = g ua, y a we have 

= -r; a i a (9) 

The non commuting operator so defined determines the parallel transport of a basis 

vector 7". It should be distinguished from the ordinary — commuting — partial derivative 

7 v jAt as defined in eq. (|37j) . Instead of the symbol <9 M Hestenes uses D M , whilst Wheeler et. 

al. JH] use V M and call it "covariant derivative" . We find it convenient to keep the symbol 

for components of the geometric operator d = 7 M <9 M . When acting on a scalar field 

2 Also the usual imaginary number i, satisfying i' 2 = — 1, is not considered as a matrix, although it can 
be represented as a matrix of real numbers. 
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the geometric derivative <9 M happens to be commuting and thus behaves as the ordinary 
partial derivative. 

Using (JSJ), eq.(0) becomes 

da = Ylu(d^ + r>°) = 7 M 7,DX = 7YD,fl, 

= 7^ • 7^D M a y + 7 M A (10) 

where 7 M • 7^ is the symmetric and 7^ A 7^ ee |(7 m 7 !/ — 7^7^) is the antisymmetric part 
of the Clifford product 7^7". In general we have 

d n a = yi 7 w... 7 ^ 7 ^ DmD/i2 ...D Atn a I/ (11) 

Here D M is the covariant derivative of the tensor calculus with the properties: 

D M = <9 M 

etc . (12) 

From the relation (fTTj) we find that after applying the vector derivative on a scalar 
twice we obtain 

dd<P = ( 7 ^ M )( 7 ^)0 = 7 M 7 V D M D^ 

= 7^7^0,0 + 17^7^0,-0^)0 (13) 

If the connection is symmetric, i.e., when torsion C p = Y p — is zero, the commutator 
of the covariant derivatives acting on a scalar vanishes, [D M , D,]0 = — C p = 0, and we 
have 

«9<90 = sTVPut = D M D^0 = -^d^^\g^d v cj>) (14) 

which is just the D'Alambert operator in curved space. Here, as usually, g = detg^ 
denotes the determinant of the metric. 

A more elaborated discussion of the geometric calculus the reader can find in refs. 
[T2*l Here let me just mention that in geometric calculus the sum of scalars, 

vectors 0^7^, bivectors a pu ^^ A 7,, etc., is as legitimate operation as is the sum of a real 
and imaginary number. 
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This formalism should not be confused with a different formalism discussed by [20] 
and [21], where not only Clifford numbers, but also differential forms take place. 

The equations of geometric calculus considered above do not rely on a particular 
representation of Clifford algebra and thus hold in space of any dimension, including 
dimension 1. In the latter case the set of basis vectors consists of only one element 70 
whose square gives a single component of the metric g 00 — 7o • 7o — 7o which in general 
depends on position. By a suitable coordinate transformation it can be transformed to 
9oo = (dx/dx') 2 g 00 = 1. 

In this paper we do not consider the action of the derivative operator (JSJ) on a spinor. 
This could be done straightforwardly by taking into account the well know fact that 
spinors can be represented [T2] 113] IT6*] 122] as members of the left or right ideals of Clifford 
algebra. 

3 Vector integral in geometric calculus 

Before going to physics, more precisely, to quantum mechanics in curved space, some 
more mathematical preparations are necessary. First let us recall that for an arbitrary 
vector with components the expression 



has no geometric sense, since it depends on a chosen coordinate system and has no definite 
properties under general coordinates transformations. 

On the other hand, if we integrate not components a? but a vector a = a fl '~f ll we obtain 
an expression which is invariant under arbitrary coordinates transformations: 



This is so, because a vector a = 0^7^ by definition does not depend on coordinates. Under 
a change of coordinates the transformation of components 



is compensated by the corresponding transformation of basis vectors 




(15) 




(16) 



dx 13 




dx'v 
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so that the vector a = 0^7^ = a' M 7^ remains invariant. Since the volume element y\g\ d n x 
also is invariant, the "average" a in eq. (jl6|) is an invariant quantity. 

A question arises at this point of whether the integration over a vector field is a 
well defined operation from the geometric point of view: the result should be a geometric 
object — a tangent vector — of the considered curved space. We shall show that this can 
indeed be the case. Namely, if we have some rule for bringing together vectors at different 
points and integrate them at a chosen point of the manifold, then such operation is 
geometrically well defined. We will see that the integral (jlfij) . if properly interpreted, can 
in fact incorporate such a rule. 

Let us now examine the integral (JTHj) more closely. First let us introduce at every point 
x a local "Lorentz" frame 3 spanned by a set of basis vectors 7 a and their duals 7" = r) ab ^fb 
satisfying 

7a • lb = Vab , 7 a ■ ^ = (19) 

where r] a b is the Minkowski metric and r] ab its inverse. At every point the basis vectors 7^ 
can be expressed in terms of j a , and vice versa: 

7„ = (7m • 7°)7a , 7a = (7a "7^ (20) 

where 

7/, • 7" = e% , 7 a • I 11 = ef (21) 

are the vielbein field and its inverse, respectively. Inserting ()2()j) into the integral (fT7)|) we 
obtain 

0= [ d n xJ\g~\a»{x)e%{x) la {x) (22) 

The above expression f!22l) is invariant with respect to arbitrary general coordinate 
transformations of spacetime coordinates x M (passive diffeomorphisms) . In addition, it is 
invariant with respect to arbitrary passive local Lorentz transformations (which act on 
the index a) , since any transformation of e a M is compensated by the corresponding trans- 
formation of 7 a so that that the integrand remains the same vector at x. Consequently, 
the integral itself remains invariant. 

3 If we interprete our formalism as describing the non relativistic theory in an n-dimensional curved 
Euclidean space, then 7° at every point x span a local frame of a flat space tangent to x. Then instead 
of a local Lorentz transformation we have a local rotation. 
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Choice of a local Lorentz frame at every point x is in fact a choice of gauge, namely 
a local Lorentz gauge. The result of the integration (j22j) does not depend on choice of 
gauge. Therefore, in order to perform calculation, we are free to take whatever gauge 
we find convenient. This will have no influence on our result. If 7 a are represented as 
matrices satisfying (|19|) the result of the integration is a matrix a, and the matrix is the 
same regardless of which local Lorentz gauge we choose in the integral (|22|) . Does the 
matrix a represents a vector of our manifold? In other words, can a be expanded in terms 
of basis vectors 7 a or 7^ taken at a chosen point x' within a domain Q of the manifold 
over which the integration is performed? 

Let us choose a gauge in which 7 a are constant 4 (see eqs. (J3~7j) - (J3T?|0 at every point x. 
Then they can be taken out of the integral and so (|22j) becomes 

a = (jf <TxJ\g~\a»(x)e\(x)} la = A a la (23) 

The result of the integration (}2~3*j) over a chosen domain Q in the manifold is the object 
a = A a 'j a . If the Lorentz basis vectors 7 a are represented as Dirac matrices in flat space, 
then also a is a matrix. Since it can be expanded in terms of a complete set of basis 
vectors, a itself is a vector. 

Although the integral (|23|) is invariant with respect to passive local Lorentz transfor- 
mations which change e a M and 7 a , there is still a freedom to perform at every point x an 
active Lorentz transformation 

e\{x) - e'%{x) = L\{x)e\{x) (24) 

whilst keeping 7„ fixed. Under the transformations (j24j) the metric does not change: 

% • l'v = 7m • lv = g^u (25) 

Here 7^ = e a ^ a and 7^ = e' a A1 7 a are two different Clifford numbers that can be 

represented by two different matrices which solve the Clifford algebra relation ()25|). The 

integral (|TK|) or its variants ((221), P^l depend on choice of 7 M (x), i.e., on choice of e a ^(x) 

and 7 a (x). Since the choice of 7 M (x) is arbitrary, it seems at first sight that there is 

no way of defining the integral (|TB|) uniquely. But is we think deeper, we observe that 

Constructing a curved space expression by means of constant fiat space matrices 7 a and vielbein e a ^ 
is a standard procedure (c.f. the action for the Dirac particle in curved spacetime). 
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in flat space the integral, if written in the general form (|22|). is not uniquely defined as 
well. In flat space we also can distinguish between the coordinate frame 7 M , satisfying 
In " lu — 9fiv( x ), an d the local Lorentz frame j a , satisfying 7 a • 7?, = r] ab . There is still a 
lot of freedom in choosing a representation for 7^ and 7 a at every point x. But there is a 
choice which distinguishes itself from other possible choices. In flat space such a choice is 
straightforward. My proposal is now that when considering the integral in curved space, 
choice of 7^ and 7 a has to be such that in the flat space limit eq. (J22J) gives the well 
known result for the integration in flat space. In the following I will discuss such a choice. 

Let us consider a curved space V n and let us choose a point x' and the quantities 
e a fJ ,(x'), J a ( x '), 7m( x ') = e ° /j,( x ')la( x ') at that point. Then at every other point x 6 Q 
we choose n vectors 7 a (x) = e a A1 (x)7 M (x), a = l,2,...,n, such that they will be parallel 
to the vectors 7° (2') = e a fJi (x')j ,J '(x') in the sense of the parallel transport along the 
geodesic 5 joining the point x and x'. Under the parallel transport the components e a fl of 
7° transform according to 

e\{x)=g»{x,x')e a v {x') (26) 
where g ll u {x,x') is the parallel propagator [2*5] . Inverting (J26|) we find 

g»{x,x') = e\{x)e a l/ {x') (27) 

Further, we represent 7 a (x') and j a (x) by the same flat space matrices so that 

7 a (x) = 7 V) (28) 

From eq. (J28|) we have 

e%(x)Y(x) = e%(x')Y(x') (29) 

which gives 

Y(x) = e a v {x)e\{x')^{x') = g\(x, x')r{x') (30) 

The latter restriction on ^"{x) is a consequence of choice ()28p. 

5 We assume here that for any two given points within the domain fl of integration there exists a 
geodesic, and that it is unique. The cases where this is not true, are excluded from our consideration. 
We assume that in a curved space which allows for non existence or non uniqueness of geodesies, the 
integration domain has to be narrowed down to such extent, that the above complication no longer occurs. 
We do not claim that our integral can be defined globally. 
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Inserting (J2fij) and (J2~H|) into the integral (J22j) we find 

= J d n x^\a^(x)g^(x,x')j v (x') 

= Jd n x^\g~\a u {x',x)^ u {x') = A v {x') lv {x') (31) 

where 

a v (x',x)=a^x)g^(x,x') (32) 

This is the relation for parallel transport of a vector with components a^(x) from a point 
£ along the geodesic to a point x' (see, e.g., fZfl). 

In eq. (|31|) vectors are first parallelly transported from points i 6 fi to a chosen point 
x' and then they are integrated. This obvious result, that also holds in flat space, has 
been obtained from the general expression (J22)) for the choice of e a ^ as given in eq. (|26j) 
and following the procedure of eqs. (|2~% j) -(|3*T ]) . 

Clearly A" (x')'y v (x') in eq. (|31|) is the vector at x' that is obtained by summing (more 
precisely, integrating) the vectors 

aV, x) lv {x') = a»{x)e\{x)e a v {x') lv {x') (33) 

which are all taken at the same point x' (which can be any point within the domain Vt). 
The argument x of the integration determines the points from which the vectors a u {x' , x) 
were brought by means of the parallel transport along the geodesic joining x and x'. The 
integral (|16j) was thus shown to have a well defined geometrical meaning. The integral 
by its construction and the choices \2b\) . \2§) automatically implies the rule for bringing 
together vectors at different points. So far it has been taken for granted that in the integral 
such as (I16|) we are summing vectors at different points and therefore such an operation 
has no geometric sense. It has not been realized that if in the definition of the integral 
(fTBj) we employ the choices (j2T?j) . (J2~%j) . then vectors are actually integrated at the same 
point of the manifold. In curved space we are doing exactly the same as in flat space: we 
first bring vectors together by parallel transport and then perform the summation or the 
integration. Eq. (jTfij) is nothing but a short notation for such operation. 

We have shown that the integral (|16|) can be written in terms of 7„ that are "constant" 
at every point x. In order to understand what a constant 7 a precisely means, let us recall 
that we distinguish between two different types of derivative: 
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(i) Geometric derivative (determining the parallel transport) gives 

dul, = r*,7, (34) 

dvla = ~^avlb (35) 

d u e% = T^e% - u ab u e b , (36) 

where u> ab ^ is the connection for the local Lorentz frame field 7 . 
(ii) The commuting partial derivative satisfies 

Tav = ^ p Vi T1 P - u abv e b ^ a + e a v ria,v (37) 

7o,i/ = arbitrary (in general may be different from zero) (38) 

,,, -V^e%-u a \e hli (39) 



Since e a ^ are the scalar components of the vector r y a = e a At 7 M , the geometric and the 
partial derivative of e a M coincide. 

Eqs. f)34j l -(|39 p hold for arbitrary 7 a and they are covariant with respect to general co- 
ordinate transformations (|17| ) .(|18 |) and local Lorentz transformations 7^ = L b ^b provided 
that the connections transform as 

IP dx' p dx a dx 13 ^ a dx' p d 2 x a 

Uf * ~ 'dx^dx^dx^ a/3 + dx a dx'»dx w ^ ' 

u' ab ^ = u cd ^L a c L b d + L a c L cb tfl (41) 

If la,u — 0, a = 1,2, ...,n, v — 1,2, ...,n, i.e., if the (commuting) partial derivative of 
7 a is zero, then we say that 7 a is constant at every point x of the manifold. 

Taking the geometric derivative of an arbitrary vector A we have after using eq.(|35j): 

d v A = d u (A a la) = {d„A a + uj\ v A b ) la = D,A a la (42) 

On the other hand, using eq. (|54*jl . we have 

d v A = d v (A^) = {d v A» + Y% A p )^ = D„A^ (43) 

where D U A P = A^. u = d v A^ + T p p A p is just the covariant derivative of a vector field. 
Comparing (|42|) and (|43|) we find that 

D u A^ = D„A a la (44) 
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Taking the commutators of the geometric derivatives acting on A = A^^y^ we obtain 
[d a ,d p ]A = A»[d a ,d p ]^ = ([D a ,D p ]A^ = R^vap A v 7^ (45) 

where 

= d^ta - doFup + ^au ~ ^ap ^ (46) 

If acting on A = A a, j a , the commutator gives 

[d a , dp\A = A a [d Q , df,) ltt = (p a , Dp]A a ) la = R aba/3 A b j a (47) 

where 

Raba/3 = — {^abafi ~ ^abf3,a + ^aca^bfi ~ ^acfi^ba) (48) 

If, in particular, we take for A just a vector j a = e^lp-, a = 1,2, ...,n, we find from 
IfEZjl that 

e/ ; „ = e a % + r^e/ = -u a \e b v (49) 

which is consistent with eq. f!36|) or (j3T?j) for the derivative of vielbein. So we have verified 
that the relation (j2Sj) for " constant" 7 a is consistent with eqs. (IM Jl - (p$fi|) for geometric 
derivative. If in a curved manifold the field j a (x) is constant in the sense that 7 QjM = 0, 
then this reflects merely a choice of the representation for 7 a (a solution to the Clifford 
algebra relation ()19|1 ) at every point x. We have shown that when expanded in terms of 
the basis vectors 7 M (x) according to 7 a = e a M 7 M the vector field j a behaves correctly from 
the geometric point of view. If vectors r y a (x), a = 1, 2, n for any x G £1 are parallel (in 
the sense of parallel transport along a geodesic) to corresponding vectors j a (x') at a given 
point x', this does not imply that vectors at different points x ^ x' are parallel amongst 
themselves. 

Extrinsic and intrinsic integral In our approach only the intrinsic geometry of the 
manifold has been considered and the integral has been given the geometric meaning 
which is intrinsic to the manifold. Had we considered our manifold V n as being embedded 
in a higher dimensional space Vn, then we would have two options of how to perform (or 
better, to define) the integration: 

(i) either intrinsically, by using the intrinsic basis vectors 7 M and 7 a satisfying the 
relations ()34 p -([35| ) and choosing a gauge in which 7 a are constant according to eq. (|35jl . 
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(ii) or extrinsically, by considering 7^ and 7 a as being induced from the corresponding 
vectors in Vn and taking the analog of the relations (|51)) - (|33|) in the embedding space 
V N . 

That is, when we write an expression such as (JT5j) we have to specify (e.g., by using 
an extra label) over which manifold the integration is to be performed. So we define two 
different integrals: 

= vector in V n (50) 

v n 



f d n x y[\g\ = f d n x J\g\ a»e\ la 



A = [ d n x = I d n x sMa^ V M j M 

J v N J 

= J d n x^\a»d,i 1 M E A MlA 



v N 

= vector in Vn (51) 

v N 

Here r] M (x^), M = 1,2,...,N, are the embedding functions for the manifold V n which 
is considered as an n- dimensional surface embedded in Vn- The quantities 7m, 1a and 
E A m = r y A ■ 7m they all refer to the embedding space Vn and are the basis vectors, the 
local Lorentz vectors and the vielbein, respectively. 

In eq.([50 )L ([51| ) l^\ v and 7 M | v are totally different objects: the former are generators 
of the Clifford algebra of V n , whilst the latter are the linear combinations of the generators 
7m of the Clifford algebra of Vn- If represented as matrices, then 7 J and 7 J are 
completely different kinds of matrices, although they both represent the tangent vectors 
to V n . This clarifies why the vector integral such as (|16|1 can be consistently defined 
either as a geometric object residing in the considered manifold V n , or alternatively as a 
geometric object residing in the embedding space Vn- 

As an example imagine a curved 3-dimensional surface V3 embedded in flat 4- 
dimensional space(time) V4. There are two possibilities: 

(i) Tangent vectors of V3 can be expanded in terms of three position dependent 2x2 
matrices a a , a = 1,2, 3, which in turn can be expressed by means of the "dreibein" in 
terms of 2 x 2 Pauli matrices <7j, i = 1, 2, 3. Summation (or integration) of those tangent 
vectors gives an object which is itself a tangent vector. 

(ii) On the other hand, tangent vectors of V3 can be induced from V4, i.e., they can be 
considered from the extrinsic point of view and expanded in term of the 4x4 matrices 
7 M , n = 0, 1, 2, 3. It is clear now that in this (extrinsic) case the summation (or the 
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integration) of the tangent vectors gives a vector in V4 which is not necessarily a tangent 
vector of V3. 

Whilst in Case (ii) tangent vectors are 4x4 matrices that can be expanded in terms 
of matrices 7^, in Case (i) there is no way to express the tangent vectors — represented as 
2x2 matrices — in terms of 4 x 4 matrices 7^. In the case (i) tangent vectors are true 
intrinsic vectors to V3 and they bear no relation to the embedding space. 

Summary - We are performing quite a legitimate operation of transferring vectors 
a(x) form x G Q into into a chosen point x' and integrating them at x'. As a result we 
obtain a vector a(x') at x'. 

By means of geometric calculus based on Clifford algebra such operation can be defined 
by employing the notation 



where 7 M (x) = e a fl (x)j a (x) are expanded in terms of vectors 7 a (x) which are all chosen 
to be parallel to 7 a (V) at x'. Mutually, of course, they are not parallel. For simplicity 
reasons we then omit the subscript {V n ,x'} at the integration symbol. 

Choice of the local Lorentz frame field j a ( x )i or ; equivalently, e a /1 (x) is arbitrary. 
Choice must be such that the integral coincides with the integral in flat space. Also in 
flat space it depends on choice of y a (x), that is, e a At (x). We choose vectors of the field 
7 a (x) all parallel to a given vector y a (x'), and obtain the usual result for the vector integral 
in flat space. The same choice of ■y a (x) (that is, of e a /J (x)) we keep in curved space. If 
our space V n is embedded in a higher dimensional space Vat, it must also be specified 
with respect to which space, V n or Vn, the parallelism of j a { x ) an d la( x ') is taken. The 
parallelism of 7 a (x) and r y a (x') with respect to V n defines the integral which is different 
from the integral in which the parallelism of r ) a {x) and 7 (#') is taken with respect to Vn- 

To sum up, the vector integral ()22|) is a functional of the frame field j a ( x ) {° T , equiv- 
alently, e a ^(x)). In other words, the integral a is an object which is defined with respect 
to a chosen local Lorentz frame field. It is convenient to choose the frame field such that 
in the transition from a curved to flat space the curved space integral coincides with the 
usual flat space integral. 




(52) 
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4 Fock-Schwinger-DeWitt proper time approach to 
quantum theory in curved space and geometric cal- 
culus 



When discussing the problem of quantum theory in a curved space V n the authors usually 
start from the following classical action 



where A is a fixed constant and r an arbitrary parameter, whilst , \x = 1, 2, n, are 
r-dependent functions denoting position of the "particle" in an n-dimensional space V n . 
Usage of the above action for description of a relativistic particle has to be taken with 
some caution. The topics has been extensively discussed in the literature. In relation 
to the usual relativistic theory which starts from the well known minimal length action 
and in which momentum is constrained to a mass shell, there are at least three different 
possible interpretations of the classical and quantized theory based on the action (J53J). 
According to one interpretation (J53)) is a gauge fixed action (i.e., an action in which 
reparametrization of r is fixed). It is a specail case (for p = 0) of the Schild action 
[21] which can be used for description of p-branes. According to another interpretation, 
(J53j) is an unconstrained action, analogous to the non relativistic action. A considerable 
number of authors has pursued such approach [T"""l "2*3 E6J and has provided the arguments 
why (I53*j) and its quantization is good for description of relativistic particles. In relation 
to the Fock-Schwinger proper time formalism [5j that was pursued by DeWitt [21 El [2], 
in eq. (|53|) we have nothing but an auxiliary, unphysical action, whose (e.g., canonical 
or path integral) quantization employs the unphysical Hilbert space and the unphysical 
"evolution" parameter r (which, in particular can be the proper time s). Physical (on 
mass shell) states are obtained by integrating the unphysical, r-dependent states, over 
r. In path integral quantization of ()53|) one obtains unphysical, r-dependent, Green's 
function from which, by integration over r, one obtains the physical Feynman Green's 
functions. 

In order to disentangle the problem of quantization in curved space(time) from the 
intricacies of the relativistic theory whose quantization has to take into account the mass 
shell constraint, many researchers have so far chosen first to tackle the easier problem 
which has roots in the unconstrained action ()53|) . The latter action was employed by 




(53) 
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DeWitt El H] i n his curved-space generalization of the Fock-Schwinger proper time 
technique for finding the Feynman Green's function G(x,x'). 

The canonical momentum belonging to the action (|53|) is p^ = dL/dX^ 1 , and its square 
p^ = M 2 is an arbitrary constant of motion. The Hamiltonian is 

H=^p»Pv = ^p 2 (54) 

where p = 7 M p M is the vector momentum. 

Upon quatization x^ and p M become operators satisfying 

[x»,p,]=t5» u (55) 

[x»,x»} = 0, [p„p„} = (56) 

Eigenvectors \x') satisfy 

x»\rf) = x'^x') (57) 
and they form a complete eigenbasis. They satisfy the following normalization condition 

urn 

(x\x>) = = ^=f2 5(x, x') (58) 

In this representation the position operator is diagonal, but the momentum operator is 
not diagonal. We shall now determine its matrix elements. 

Multiplying from the left by basis vectors 7^ and summing over v we obtain the 
commutator which involves the vector momentum p: 

[x",p]=i'f (59) 

If we sandwich the commutation relations ()59j) between the position eigenstates we 
have 

(x\[x^,Y(x)Pu]\x'} = (x» -x'»)(x\YPu\x) = (x\iY\x'} = iY(x)6(x,x') (60) 

From 



j a {x)d a ^{x fl 5{x,x')) = 1 a {x)^ a 5{x,x') +x^d a 5{x,x') 
1 a {x)d a [x' IM 5{x,x')) = -i a (x)x'^d a 5(x,x') (61) 
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we obtain 

(a? - x'^ a (x)d a 5(x, x') = -^5(x, x') (62) 
Comparing eqs.(JB0J) and (jHSJ) we obtain 

(x\YPuW) = -ij v dj{x, x') + F(x)8{x, x') (63) 

where F(x) is an arbitrary position dependent vector. 

The same relation ()63|) can be derived from the commutation relations (fo3j) which give 

(x\p^\x') = -id^6(x, xTj + F^(x)S{x, x') (64) 

and contracting eq. (JMj) by 7 M . 

In the literature so far the authors considered the matrix elements (I64|) and required 
that they had to satisfy the condition for Hermiticity 

(x'lp^x)* = (x\p^\x') (65) 

This together with the commutation relations [p^,^] = restricts the choice of F^x): 

F M = fy(~ln|0|-x) (66) 

where xi x ) is an arbitrary function - a phase. Hence, for the choice x = 0> 

(xlp^x'} = -i (dp + jdfjln \g\J 8(x, x') (67) 

where 

Rln \g\ = \\g\-%\g\ = Wg^d^ 12 = l^Wo^ 4 = (68) 

By employing the powerful geometric calculus we recognize that the momentum op- 
erator should be defined as the vector operator 

p = 7 % (69) 

So instead of the matrix elements (xlp^x') (eq. ()64)l ) we have to consider the matrix ele- 
ments (x\p\x') (eq. ()63)l ). Hermiticity condition for the latter matrix elements determines 
F(x) in eq-dSU). It turns out (see eas. (f73jl - (f7Hj) l that 

F(x) = (70) 
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Hence 

(x\p\x f ) = -i^(x)d„5(x,x') (71) 

and we have 

(x\p\<f>) = J (x\p\x')^\g{x')\dx' (x'\<t>) = -i^d^ = p<\> (72) 

where \(f>) is a state vector and (x\4>) = a wave function. 

We shall now demonstrate that the matrix elements (j71j) are Hermitian. Using 

d,5(x, x') = -& 6{x, x') - -=L= d'^^)\5{x, x>) (73) 

which follows directly from eq. (JEHJ) , and by taking into account eq. (}8T|) (see next section) 
and the relation 

1%, = -i!=r wsi (74) 



we have 

- iY(x)d fl 5(x, x') = i^(x')d'^(x, x') (75) 
From the latter relation it is straightforward to verify that 

{x'\p\x\)* = (x\p\x') (76) 

Thus the matrix elements ()71|) of the vector momentum operator in curved space satisfy 
the Hermiticity condition. No extra terms to eq. ()71|) are necessary in order to assure 
Hermiticity. This is not the case in the approaches which work with components p^ of 
momentum operator: the matrix elements (xlp^x') = —id^5(x,x') are not Hermitian. 

Let us now also calculate the matrix elements of the square of the vector momentum 
operator. Using eqs. (f7T|) . (f75]) and (J8§|) . namely d^^Jlg]^) = 0, we find 

(x|p 2 |0) = J (x\p\x') <J\g{x')\d n x' (x'\p\x") ^f\g(x")\d n x" (x"\(f)) 



xiY(x")d"( 5 ^' X "n J\g{x")\d n x"<p{r,x" 

' 1 WOI / 



i [ 7*00% [%==] v^(^ dV 
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where <f>(r,x) is the projection of a state \<p) onto position eigenstates. The above result 
explicitly demonstrates the consistency of our procedure within the matrix formalism of 
quantum mechanics in curved space. We see that the matrix element {x\p 2 \(p) can be 
calculated by inserting twice the complete set of the position eigenstates so that under 
the integration we have the product of two terms such as (x\p\x'). The result is the 
product of two vector differential operators —ij^d/j,, acting on the wave function, which, 
according to eqs.(J7HHj), is equal to the covariant d'Alambert operator (multiplied by —1). 

The quantization procedure introduced in this section does not use at all the vector 
integral which was discussed in previous section. Because of the presence of the delta 
function, the integrals such as ()72|) and (|77jl are local. The former integral gives a vector, 
whilst the latter integral gives a scalar at a point x. The delicate vector integral will 
be used in next section only in the calculation of the expectation value for the vector 
momentum operator. 

5 On the Schrodinger representation for the momen- 
tum vector differential operator and its expecta- 
tion value 

We have seen (see eqs. (|59^) - (|72*j) ) that in the Schrodinger (coordinate) representation 
the momentum operator, that remains Hermitian in curved space, is a vector differential 
operator 6 

p=-id = -vfd^ (78) 

6 As = —id^, also the vector operator p = j^p^ is generator of translations, since 6(f) = d^ipSx^ = 
ip^4>8x^ — i(p ■ 5x)4>, where 5x = 8x^ 7^. 
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acting on a wave function <f)(r, x^) which is taken here to be a scalar. The wave function 
is normalized according to 



\g\ d n x (j)*(r, ^)0(r, a^) = 1 (79) 
and satisfies the Schrodinger equation 



i-^ = (80) 



#0 = ^ = ^H) 2 «9 2 = -^D M D> = ~ -]L dM\g\ sT d v <f>) (81) 



where 



\9\ 

With the definition (I78j) o/ momentum operator there is no ordering ambiguity in the 
expressions such as p 2 , p 3 , etc. . We have 

p4> = -iY l d^4> 

p ^ = {-i-fd»){-ii'd v )<j>= (-OVt^d^ 

p 3 = (_i 7 ^ M )(_i 7 ^)(_i 7 ^ a )0 = (_i)3 7 ^ T a D ^ Da (82) 
etc. (83) 

All those expressions are covariant with respect to arbitrary coordinate transformations 
in curved space. The metric g^ v of curved space is implicit in the position dependent 
basis vectors 7^(2) satisfying eq.(fTJ). There is no ambiguity of where to place 7^, 7", etc., 
in the product of any number of operators p. 

Remember that we are considering here the action of the geometric derivative <9 M on 
a scalar so that <9 M coincides with the partial derivative of 0. Had we considered the 
action of on a spinor, the situation would be different. Consideration of spinors is 
beyond the scope of this paper. 

Let us now explicitly verify that the operator p is self-adjoint with respect to the scalar 
product 

r d n xJ\g\(f)*p<f) = (p) (84) 



where = 0(r, x) depends on the evolution parameter r and coordinates x M . 

Eq.(jSU) defines the expectation value of the momentum operator p. By using eq.(f78|) 
we write it explicitly 

(p) = - % \ VbK^v^ ( 85 ) 
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Its complex conjugate value is 7 

(P)* = if V^1d n x07^ t 0* 

= -if \f\g\d n x frfd^-i f d n x<f>*d li (y/\g\'f) 

+ i f d n xd^*J\g~\r<P) (86) 

From eq. we have 

= -r£ a 7 Q (87) 
On the other hand, the derivative of the determinant gives 

d,^\g\ = ^\g\K, (88) 

Therefore, 

^(VW') = >f + \f\g\d^ = ^\(T^Y - V) = (89) 

Using eq.(JEHJ) and assuming that 4>(t,x) vanishes at the boundary of the integration do- 
main Q, so that the boundary term can be omitted, we find from (JSfiJl that the expectation 
value is real: 

(p>* = (P) (90) 

which means that the vector momentum operator p = — ry M <9 M is self-adjoint with respect 
to the scalar product 

The boundary term in eq. ()86|) comes from the boundary theorem which for an arbitrary 
polyvector A in curved space V n can be written as |T2| 

/ du-3A = ( da A (91) 

JQ JT, 

Here Q is a closed r-dimensional volume bounded by an (r — l)-dimensional surface S, 
whilst the r- vector du is the volume element of Q, and the (r — l)-vector da is the surface 
element of E. 

It is straightforward to show that in Riemanian geometry the boundary operator d 

in eq. (J91|) may be taken to coincide with the geometric derivative defined in eqs. (JSjHQ 

7 The complex conjugate vector a* is defined according to the relation a* = a* M 7 M in which the basis 
vectors remain unchanged. 
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and used in eqs. (|85j ) -([90 p . It cannot coincide with the commuting partial derivative with 
the properties flHUHflUHJ). 

The expectation value (|85|). which we will now denote (p) = P, is defined by the vector 
integral which includes into its definition a choice of a point x'. Since the wave function 
depends on time r, the expectation value (p) = -P(r) in principle also depends on r. 

Let us now calculate the derivative of the momentum expectation value with respect 
to r. From 

(p) = J <j)*p(f)^\d n x (92) 

we have 

d(p> 



dr 



/ (^W + <^^) = * / (WW - <f>*H(f>) ^\d n x (93) 



where we have taken into account the Schrodinger equation idtfi/dr = H<p and its con- 
jugate —id(j}*/dr = Hip*. It is straightforward to show that for H (which is Hermitian, 
since p is Hermitian) we obtain 

J (H<f>*)p<f> ^\d n x = J <f)*Hp<f) sj\^\ d n x (94) 



so that eq.(jH3J) becomes 

d(p) 



[H,p}0J\g\d n x = O (95) 



dr 

due to the fact that p commutes with H = Ap 2 /2. So we have found that the expectation 
value (p) does not change with r. 

We can also look at the r derivative of (p) from another angle. According to our 
interpretation of the vector integral as discussed before (see eqs. (ITB jl - ()l*Sj )). the result of 
the integration is a vector at a chosen point x' of the manifold. In the following we will 
omit the prime and denote the chosen point by x. A question arises as to which point 
x we should choose. A natural choice for x is the expectation value (x) = X(t) of the 
particle's position. The latter point changes with r, therefore at every r the momentum 
P(t) is taken at different point X(t). We can compare P at different points by means of 
the geometric derivative which performs the parallel transport of P at x^ + 5x^ to the 
point x M so that 

P(x + 8x) - P(x) = 5P = d^P 5x" (96) 
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By taking x^ = X M (r) and x^ + Sx^ = X M (r + St) we have that Sx^ = X^St so that 

SP = O i1 PX»St (97) 

In (19 7j) we calculated the change of P along the trajectory X m (t) at two infinitesimally 
separated times r and r + St. The derivative of P is thus 

dP _ d(p) 



dr dr 
Using and we thus have 



^PI" (9* 



dP^ 



x^p = x^{P v lv ) = x^p* + r^p")T, = \ -^r + r; p x^p" j T „ = o (99) 

By putting P u = (1/A) dX^/dr, where X u = (x u ) is the expectation value of the particle's 
position (i.e., the "center" of the wave packet), we find that (|9*9*j) is the equation of motion 
that follows from the classical action f)53jl . Since, according to (|95jl the momentum vector 
P is a constant of motion, also its square P 2 = P V P V = M 2 is constant. This gives the 
relation A~ 2 X^X M = M 2 so that P v = M X v /(X^X^) 1 ' 2 . Inserting the latter expression 
for P v into eq. ()99)1 we have 

i <• ' *" \ r = o (1M ) 



v^2 dr Vv 7 ^ 2 "/ X 2 
which is i/ie geodesic equation. So we have found a natural result that the expectation 
value of the momentum operator follows a geodetic trajectory in our curved space. That 
is, at every r, the expectation value P is tangent to a given geodesic. 

6 Conclusion 

We have shown how the long standing problem of the ordering ambiguity in the definition 
of the Hamilton operator for a quantum point particle in curved space can be elegantly 
resolved by the quantization which employs the geometric calculus based on Clifford 
algebra. This is yet another besides at least two other approaches to quantization in curved 
spaces, namely "geometric quantization" [II. and Kleinert's path integral quantization 
[7j, in which ordering ambiguities do not arise. In most other formulatations of quantum 
mechanics in curved spaces ambiguities occur due to arbitrary choice of operator ordering 
prescription. Common for those formulations is that they do not use a geometric language. 
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Instead, they use the component, coordinate-based, notation of tensor calculus, which is 
a very powerful mathematical tool, but also has its limitations. Our conclusion is that 
momentum and Hamilton operator in curved space cannot be consistently formulated 
in terms of the component notation. The geometric calculus based on Clifford algebra 
provides a natural definition of the momentum operator p which is Hermitian, and of 
the Hamilton operator which is free of the ordering ambiguities. Moreover, the language 
of geometric calculus enables us to define and calculate the expectation value of the 
momentum operator which turns out to follow a classical geodetic line — a very reasonable 
result. We have also been able to handle consistently the matrix elements of p and p 2 
between the position eigenstates, by which we have further completed the formulation of 
quantum mechanics in curved space. This opens new perspective on the subject and its 
further development. 
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